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ABSTRACT 



Aims. We aim to study the structure of the nebula around the post-AGB, binary star 89 Her. The presence of a rotating disk around 
this star had been proposed but not been yet confirmed by observations. 

Methods. We present high-resolution PdBI maps of CO J=2-l and 1-0. Properties of the nebula are directly derived from the data 
and model fitting. We also present N-band interferometric data on the extent of the hot dust emission, obtained with the VLTI. 
Results. Two nebular components are found: (a) an extended hour-glass-like structure, with expansion velocities of ~ 7 km s -1 and a 
total mass ~ 3 10~ 3 Mq, and (b) an unresolved very compact component, smaller than ~ 074 and with a low total velocity dispersion 
of ~ 5 kms~'. We cannot determine the velocity field in the compact component, but we argue that it can hardly be in expansion, 
since this would require too recent and too sudden an ejection of mass. On the other hand, assuming that this component is a keplerian 
disk, we derive disk properties that are compatible with expectations for such a structure; in particular, the size of the rotating gas disk 
should be very similar to the extent of the hot dust component from our VLTI data. Assuming that the equator of the extended nebula 
coincides with the binary orbital plane, we provide new results on the companion star mass and orbit. 
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1. Introduction 

Planetary and protoplanetary nebulae (PNe, PPNe) very often 
show axisymmetric shapes and fast axial expansion, which are 
thought to be due to shock interaction between the very col- 
limated post-AGB jets and the slow and isotropic AGB wind. 
Theoretical calculations show that accretion from rotating disks 
onto the post-AGB star or a companion can explain the axial 
jets (see lSokerfl2002tlFrank & Blackmanll2004 etc) and provide 
the high energy and momentu m required to explain the PPN dy- 
namics (Buiarrabal et al. 2001). On the other hand, the presence 
of a large amount of orbiting material has been invoked to ex- 
plain some remarkable observational properties: 1) the peculiar 
abundances sometimes found in the atmospheres of post-AGB 
stars, where refractory elements are deficient, presumably be- 
cause of reaccretion of material in which grains have been effi- 
ciently form^djremoving such elements from the gas) and ex- 
pelled dVan Wincke 11120031) ; and 2) the detection of a NIR excess 
due to dust at temperatures close to that of sublimation, which 
should correspond to stable structures (i.e. in keplerian rotation) 
close to the star. Such SEDs are r ather common among post- 
AGB stars and, in a review paper of de Ruv ter et al.l d2006), fifty 
one objects are listed. These features are systematically found 
to be associated to binarity, which would explain the excess of 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer, as well as on observations of the Belgian Guaranteed 
time on VISA (ESO). IRAM is supported by INSU/CNRS (France), 
MPG (Germany) and IGN (Spain). 



angular momentum needed to form keplerian disks from previ- 
ously ejected gas. 

The existence of rotating disks around post-AGB stars is 
therefore a basic question to understand the post-AGB evolution, 
in particular, the shape and dynamics of PPNe. 

Disks or tori of molecular gas around post-AGB stars are 
commonly detected as the central part of PPNe, although in gen- 
eral they are not observed to rotate, but to be in exp ansion, like 
the re st of the nebula. It has been suggested (e.g. Fon g et all 
12006b that wide oblate structures found around a few AGB 
stars could be in rotation, but the evide nce is very controversial 
dJosselin et alJl2000l:lHirano et al.ll2004 etc). Direct detection of 
a keplerian velocity field has been obtained for o nly one PPN till 
now: the Red Rectangle (Buiarra bal et al.l l2005). a well-known 
PPN surrounding a binary star. 

89 Her (HR6685) was first studied by iBTdelmanl dl95lh . who 
noted the presence of a presumably massive supergiant at high 
Galactic latitude. There is now, however, general consensus that 
89 Her is a low-mass object in its post-AGB phase of evolution. 
It is one of the few post-AGB stars with a significant Hipparcos 
parallax, placing it at 1^ kpc from the Sun. Interestingly, 89 Her 
is a binary syst em that shows both a conspicuous NIR excess, 
due to hot dust (I Waters et al.ll 1 99 3b Ide Ruvter ei~al"1l2006l). and a 
far-IR flux indicating the presence of large grains ( Shenton et al. 
1995). Its CO lines show a peculiarly sh arp profile very simila r 
to the one detected in the Red Rectangle (Buiarra bal et al.l2001l) . 
Previous interferometric maps of the CO emi ssion from 89 Her 
(lAlcolea & Bujarrabal|[l99"5l;lFong et alj|2006l) have not yielded 
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Fig. 1. Channel maps of the 12 CO 7=2-1 (upper panels) and /=l-0 lines (lower panels) from 89 Her. The first contours are, 
respectively, 15 and 30 mJy; the contours are separated by a factor 2 and the negative contours (at -15 and -30 mJy, respectively) 
are indicated by dashed lines. The LSR velocity in kms" 1 for the center of each channel is indicated in the upper left corner. The 
J2000 coordinates for the reference position, the cross in the maps, are R.A. 15:55:25.19 and Dec. +26:03:00.0. The black ellipses 
in the last panels indicate the beam half-intensity sizes. 



conclusive results, due to the small size and low intensity of the 
source. A roughly rounded structure was found, showing a cen- 
tral core plus an extended clumpy halo, with no clear velocity 
pattern. 

2. Observations 

We report on high-resolution interferometry with the IRAM 
Plateau de Bure Interferometer of the 12 CO(1-0) and 12 CO(2-l) 
lines and continuum emission from 89 Her. Observations were 
carried out under excellent atmospheric conditions (pwv = 0.5- 
2 mm, seeing 0'.'l-0'.'4) with the six-element array in the D 
configuration (November 11, 2005) and in the newly extended 
A configuration (January 19 and 25, 2006). The cross-correlator 
was set up to cover the lines with 20 and 40 MHz bands each and 
the continuum emission with two 320 MHz bands. The flux cal- 
ibration was referred to MWC349 and 3C273, while phase and 
amplitude calibrations were done using 1751+288. 

Data were calibrated in the antenna-based manner using the 
GILDASQ software package. Velocity maps were produced with 
natural weights at 2.6 mm in order to analyze the extended 
molecular distribution and by applying an inverse Gaussian ta- 
per at 1.3 mm to analyze in detail the compact emission in 
the immediate surroundings of 89 Her. Cleaned and primary 
beam-corrected maps were obtained with synthesized beams of 
3'.'60x3'.'32 at position angle -97° and 0768x0760 at -97°, re- 
spectively at 2.6 mm and 1.3 mm. 



1 The GILDAS software package is available for download on the 
IRAM web pages. 



The line emission was resampled to an effective velocity res- 
olution of 1 kms" 1 yielding a one-sigma sensitivity level of 
~ 9.3mJy/beam for the 12 CO(1-0) and of ~ 15.3mJy/beam 
for the 12 CO(2-l) line. The continuum sensitivities were esti- 
mated from the line-free emission regions of the bands to be 
~ 0.3mJy/beam at 2.6 mm and ~ 0.5mJy/beam at 1.3 mm. 
Continuum emission was well-detected at 2.6 mm and 1.3 mm 
but remained unresolved, yielding flux densities of 2.7 mJy and 
9.2 mJy, respectively, compatible with thermal dust emission. 
The size of the continuum source is smaller than 073. 



2.1. VLTIdata 

We also performed N-band interferometric observations using 
the two-beam VLTI/MIDI instrument. Three set-ups were used 
in April 2006 during Belgian guaranteed time on the auxil- 
iary telescopes (VISA). For one of the observations, the base- 
line of 16 m (projected angle = 78 degrees) was used, while 
the other two were conducted on a larger baseline of 32 m 
with similar projected angles (80, 77 degrees). The visibili- 
ties were calibrated using observations of unresolved calibra- 
tion stars (HD 168454 and HD 123139), which were observed in 
strict concatenation with the science target. These observations 
were performed in the context of a large program to systemati- 
cally study the circumstellar discs around binary post-AGB stars. 
With the 1.8m auxiliary telescopes, the sensitivity of MIDI al- 
lowed fringe detection on a correlated flux larger than 20 Jy at 
12 ^m. For the detail s on the data reduction, we refer the reader 
to lDeroo e t al. (200l). 
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3. Nebula structure: a rotating disk around 89 Her? 

As seen in our maps (Fig. [TJ, the nebula around 89 Her consists 
of two very different features, a quite extended one and a very 
compact central clump. The shape of the extended component 
strongly suggests an expanding hour-glass-like structure, whose 
axis is slightly inclined with respect to the line of sight. The 
extended component must have a maximum expansion velocity 
~ 6 - 7 km s _1 (the velocity probably increases with the distance 
to the star, in view of the variations in the image size with the 
velocity and as is confirmed by the model calculations explained 
below). For a distance of 1 kpc (Sect. 1), the typical total size of 
this extended component is estimated to be ~ 1.5 10 17 cm. 

We do not detect the extent of the compact component, 
which must be smaller than ~ (X'4, equivalent to a diameter 
smaller than ~ 5 10 15 cm. This size limit is comparable to that of 
the dust emitting at 1mm wavelength, and probably both emis- 
sions come from the same component. Note that interferomet- 
ric maps of the optically-thin dust emission at this wavelength 
strongly select the massive compact components. The total ve- 
locity dispersion is ~ 5 kms -1 , with no significant velocity pat- 
tern. A slight velocity gradient in the north-east south-west di- 
rection could be seen in our 7=2-1 maps of the central clump, 
which could be consistent with the trend found for the extended 
component. But note that contamination by the emission of the 
wide component may be relevant in this respect, as it happens 
in the 7=1-0 maps, where the extended component emission al- 
ways dominates the brightness distribution due to the relatively 
wide beam. An inn er E-W elongatio n also appears in the low- 
resolution maps by Fon g et all d2006l) . but our data clearly show 
that it is due to emission from the hour-glass-like structure. 

We have tried to reproduce our data by synthetic maps from 
an axially symmetric model nebula. The model description and 
assumptions can be found e lsewhere (e.g. Bujarrabal et al. 1998; 
Castro-Carrizo et a il 120021) . We considered an hour-glass-like 
extended structure, with very thin walls, plus an equatorial com- 
pact disk. 

In Fig. [2] we represent the model nebula shape that repro- 
duces our data, as well as the density and velocity distributions 
in the extended component. The symmetry axis of this compo- 
nent forms an angle of ~ 15° with the line of sight, i.e. 75 degrees 
with respect to the plane of the sky. In the extended component, 
the total density is assumed to be inversely proportional to the 
distance to the star, r, with n(10 16 cm) = 3 10 3 cm 4 . The ro- 
tational temperature of the extended nebula decreases linearly 
with r, between 1 1 K and 5 K. The velocity in this component 
is radial and its modulus increases linearly with r, between 2.5 
kms 1 in the center up to 8 kms -1 . The CO relative abundance 
is assumed to be constant and equal to 3 10" 4 in all the nebula. 
Note that we have tried to keep the model very simple, in view 
of the lack of observational information. 

The values of most parameters of the central unresolved 
component are uncertain. The velocity field of the compact disk 
cannot be reliably deduced from the data, as mentioned above, 
so only a velocity dispersion can be estimated for this compo- 
nent. We have taken conditions for which the 12 CO lines are opti- 
cally thick, in view of the lo w contrast shown in the line core for 
12 CO and 13 CO transitions dBuiarrabal et all 1200 ll) . Therefore, 
12 CO (and even 13 CO) lines could be opaque in it, and the de- 
rived density and mass in the central clump are only lower limits. 
Assuming high optical depth and the size limit given above, we 
can derive a kinetic temperature limit from the measured bright- 
ness, Tic > 60 K. 
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Fig. 2. Geometry and distribution of density and velocity in the 
model nebula. Note that we only obtained limits for the size and 
physical conditions of the central, compact component, repre- 
sented by the black circle. 

The total masses of both compone nts are compatible with 
those derived in lBuiarrabal et al.l d2001) from 13 CO lines (~ 10~ 2 
Mq for the central component and ~ 3 10~ 3 Mq for the extended 
one; note that we must correct the mass values in that paper for 
the different assumed distance). 

The best fitting of the 7=2-1 maps from our model is found 
in Fig. [3] Note that the model necessarily predicts quite sym- 
metric results for relatively positive and negative velocities, with 
respect to the equator, which is not exactly the case in the ob- 
servations. Therefore, our fitting can only represent a kind of 
average of the maps around the central velocities. The observed 
asymmetry could correspond to a small deviation of the symme- 
try axis for both cups with respect to the line of sight, by about 
10 degrees; a similar deviation for the projection of the axis on 
the plane of the sky is also directly suggested by our maps. 

From the size and velocity of the extended component given 
above, we find a typical lifetime ~ 3500 yr. This value is high 
compared to the typical ages of other PPNe, often around 1000 
yr, but compatibl e with expectations for post-AGB stars wi th low 
initial mass (e.g. lVassiliadis & Woodll994lBlockerill995l) . 

If we also assume that the compact component is in expan- 
sion, we can calculate a typical lifetime. If, as often happens in 
bipolar nebulae, it is a disk perpendicular to the symmetry axis 
of the extended hour-glass structure (then forming an angle of 
about 15° with respect to the plane of the sky), the lifetime must 
be shorter than ~ 100 yr. A less strong upper limit is obtained if 
we assume that the central structure is spherical and in isotropic 
expansion, in this case the lifetime must be shorter than 300 yr. 

Note the significant difference between the lifetimes of both 
extended and compact components. The very short lifetime de- 
rived for the compact core is surprising. We expect copious mass 
losses in AGB stars, but not that an F2-type star ejects about 
10~ 2 Mq in such a sudden event, completely independent of the 
ejection of the extended component. Moreover, we note that the 
star is slowly evolving now. It is known to have its present ef- 
ficient tempera ture since at least 1950 (see IWaters et alJ l 1993; 
Bidel mar3ll951h up to 2005 (unpublished data), so we do not 
expect that it was in the AGB only ~ 100 yr ago. 

The very inner nebular structure was also probed by our 
VLTI observations. The total N-band flux is completely dom- 
inated by the circumstellar hot dust component. Even on our 
longest projected baseline (31.1 m), this N-band emission is 
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barely resolved, with normalized visibility amplitudes >0.7 in 
the whole spectral band. The hot dust therefore must come from 
an angular size smaller than ~ 31 mas, i.e. a diameter 5 31 AU 
at the distance of 89 Her. This re sult is compatible wi th the un- 
resolved N-band direct image bv lMeixner et al. I (fl999h . 

If we assume that the central component is a rotating disk, 
we cannot derive lifetimes, but we can compare the velocities 
with the keplerian ones. Taking a typical projected velocity of ~ 
2 km s _1 , the characteristic rotation velocity would be ~ 8 km s _1 
for the deduced axis inclination. For a central total mass of 1 Mq 
(see mass estimate in Sect. 4), the rotating gas must be at a typi- 
cal distance of about 15 AU, yielding a typical disk size ~ 4 10 14 
cm or 30 milli-arcseconds. We note that the outer diameter of the 
rotating disk could be larger than these values, corresponding to 
lower rotation velocities but included within the observed veloc- 
ity dispersion. This very compact size is strikingly compatible 
with the size measured from our VLT interferometry at 10 pm. 
We finally note that a typical temperature for the central disk Tk 
~ 600 K is derived from the measured brightness temperature (at 
relatively high velocities but within the central component) and 
a size ~ 0705. This value is comparable to the temperature of 
the warm dust derived from analysis of the SED in 89 Her and 
to that of the keplerian disk in the Red Rectang le derived from 
CO e mission modelling, ~ 500 K at 4 10 14 cm (IBuiarraba l et al.l 
2005). 



4. Discussion and conclusions 

Our CO data of 89 Her show the presence of a double struc- 
ture, with an unresolved compact circumstellar component and 
an hour-glass-like extended outflow. Although the velocity field 
in the inner region is not well-measured, we conclude that it is 
probably a keplerian disk with a very small extent (diameter 5 
10 15 cm, £ 07 1). This is supported by our N-band interferometric 
data that shows that, in the N-band, the extent of the dust emis- 
sion is barely resolved at a 31.1m baseline. The dust emitting in 
the 1mm continuum in 89 Her is also confined to a small central 
region (diameter < 073), probably the same compact component. 

The existence of rotating inner disks has been proposed to 
explain some properties of hot dust emission from other binary 
post-AGB stars (Sect. 1), properties that are in fact shared by 
the Red Rectangle and 89 Her. Our results support the identifi- 
cation of inner stable reservoirs of material in keplerian rotation 
in these objects. The small size of the keplerian disk in 89 Her 
can be a common property of such candidate disks, in view of 
the high dust temperature measured in them. 89 Her and the Red 



Recta ngle are likely the nearest of these objects (Ide Ruvter et al.l 
2006), so that these small regions are not expected to be resolved 
by present day radio interferometers, and a direct confirmation of 
the keplerian dynamics in the other objects will require the use of 
future instruments. Our N-band interferometric experiments do 
confirm the very compact nature of the w arm dust emission in 
similar objects like HD 5298 1 and SX Cen dDeroo et al.ll2006l) . 

Assuming that the symmetry axis of the resolved hour-glass 
outflow is perpendicular t o the inner orbit of the binary, the mass 
function of 0.000838 M (IWaters et al.ll 19931) converts to a mass 
of the companion of only 0.35 M Q (i=10° gives 0.62 M Q , i=20° 
yields 0.25 M Q ) assuming a primary of ~ 0.6 M . Since there 
is no evidence of a hot, compact degenerate component in the 
system, the companion is likely to be a very low-mass main se- 
quence dwarf. 

With this inclination the orbit has a semi-major axis of 0.31 
AU (a sin i = 0.08 AU). This orbit is well within the sublimation 
radius of the primary, so the circumstellar dusty disk must be 
circumbinary. The primary does not fill its Roche-lobe now; but 
assuming the star had a similar luminosity when on the AGB, it 
must have had a phase of very strong binary interaction. 
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